%% Implied Transition Probabilities
function f= TransProb(j,E_l,params)              
    a2 = params(j,1);
    b2 = params(j,2);
    a3 = params(j,3);
    b3 = params(j,4);
    a4 = params(j,5);
    b4 = params(j,6);
    a5 = params(j,7);
    b5 = params(j,8);

    D(1) = 1;
    D(2) = exp( a2+b2*E_l );
    D(3) = exp( a3+b3*E_l );
    D(4) = exp( a4+b4*E_l );
    D(5) = exp( a5+b5*E_l );
f   = D(:)./sum(D(:));  
